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Abstract 

We determine and study the ground states of a focusing Schrodinger equation in dimension one 
with a power nonlinearity and a strong inhomogeneity represented by a singular point 

perturbation, the so-called (attractive) 5' interaction, located at the origin. 

The time-dependent problem turns out to be globally well posed in the subcritical regime, and 
locally well posed in the supercritical and critical regime in the appropriate energy space. The 
set of the (nonlinear) ground states is completely determined. For any value of the nonlinearity 
power, it exhibits a symmetry breaking bifurcation structure as a function of the frequency (i.e., 
the nonlinear eigenvalue) lo. More precisely, there exists a critical value u* of the nonlinear 
eigenvalue w, such that: if wo < < j then there is a single ground state and it is an odd 
function; if a; > w* then there exist two non-symmetric ground states. 

We prove that before bifurcation (i.e., for lo < uj*) and for any subcritical power, every ground 
state is orbitally stable. After bifurcation (a; = w* + 0), ground states are stable if fx does not 
exceed a value fi* that lies between 2 and 2.5, and become unstable for /i > /i*. Finally, for /i > 2 
and w ^ w*, all ground states are unstable. The branch of odd ground states for w < w* can be 
continued at any u) > u* , obtaining a family of orbitally unstable stationary states. 
Existence of ground states is proved by variational techniques, and the stability properties of 
stationary states are investigated by means of the Grillakis-Shatah-Strauss framework, where 
some non standard techniques have to be used to establish the needed properties of linearization 
operators. 



1 Introduction 

The present paper is devoted to the analysis of existence and stability of the ground states of a 
nonlinear Schrodinger equation with a point defect in dimension one. The Schrodinger equation 
bears an attractive power nonlinearity, and the defect is described by a particular point interaction 



1 



in dimension one, the so-called attractive 6' interaction. In a formal way, the time dependent equation 
to be studied is given by 

r idtm = - d.xm - i^'^m - \\m?^m ^^ 
I m = V'o ^ ^ 

where Vo represents the initial data, and A > 0, 7 > 0, /i > 0. The non rigorous character of the 
expression (jl.ip is due to the fact that the combination ^dxxi'it) + jd'Qip{t) is meaningless if literally 
interpreted as an operator sum or as a form sum, due to the exceedingly singular character of the 
perturbation given by 6q. Nevertheless, it is well known ([5]) that it is possible to define a singular 
perturbation of the one-dimensional laplacian — ^ which is a self-adjoint operator in L^(M) and 
has the expected properties of the stated formal expression. The self-adjointness is implemented 
through the singular boundary condition which defines the domain of the operator, i.e 

D{H^) = {i;e H\R\{0}), \ i;'{0+) = V^'(O-) , V'(0+) - V(O-) = -7^'(0+) } 

and the action is H^ip = —tp" out of the origin. 

The 5' is called repulsive when 7 < and attractive when 7 > 0. Note that the generic element of 
the domain of 5' interactions is not continuous at the origin (but its left and right derivative exist and 
coincide), at variance with the milder and better known case of the 5 interaction (which is defined 
by £ -ff"'^(M) and boundary condition '0'(O'*") — ■(/''(O") = aV(0)i where the real parameter a is 
interpreted as the strength of the interaction). The quadratic form computed on a domain element 
ip turns out to be 

which justifies the name of 5' interaction given to H^. Concerning the physical meaning of the 
5' interaction, perhaps its best known use is in the analysis of Wannier-Stark effect (see [6]). It 
is known that the spectrum of Wannier-Stark hamiltonians in the presence of a periodic array of 
5' interactions shows a remarkably different behaviour with respect to the case of regular periodic 
potential; in particular, the spectrum has no absolutely continuous part and it is typically pure point 
(by the way, the corresponding properties of Wannier-Stark for the 5 array are not known). In the 
present paper, the core interpretation is that of a strongly singular and non trivial scatterer. It is 
known that the 8' interaction cannot be obtained as the limit of Schrodinger operators in which 
the potential is a derivative of a (5-like regular function, as the name could erroneously suggest 
(see [111 [T4| for a thorough analysis of this problem). An approximation through three scaled 5 
potentials exists, but the scaling is nonlinear as the distance of the centres vanishes. Nevertheless, 
the 8' interaction has a high energy scattering behaviour that can be reproduced, up to a phase 
factor, through scaling limits of scatterers with internal structure, the so-called spiked-onion graphs 
(see [6]). These are obtained joining two halflines by N edges of length L and letting L ^ 0, N ^ 00 
while keeping the product NL fixed. An analogous behaviour is obtained considering a sphere with 
two halflines attached. These results enforce the interpretation of the 5' interaction as an effective 
model of a scatterer with non elementary structure. By the way, in this respect the results of the 
present paper give support to this view through the analysis of the bifurcation of nonlinear bound 
states. 

From a more abstract point of view, both 6 and 6' interactions are members of a 4-parameter 
family of self-adjoint perturbations of the one dimensional laplacian, the so-called 1-dimensional 
point interactions (see [HEKTS]). As explained above, we interpret the presence of a point interaction 
in the equation (II. Ih as a model of strongly singular interaction between nonlinear waves and an 
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inhomogeneity. When the inhomogeneity is described hy a S interaction, a fairly extended hterature 
exists of both physical and numerical character; more recently, there has been a growing interest 
in this model from the mathematical side, in the attempt to establish rigorous results concerning 
the existence of stationary states ( [9l [Mj ) , the asymptotic behaviour in time ( [25] ) , and the reduced 
dynamics on the stable soliton manifold ([2D]). To the knowledge of the authors the only rigorous 
result about the NLS with 6' interaction concerns the well-posedness of the dynamics, and is contained 
in [1] , where the whole family of point perturbation is treated in the presence of a cubic nonlinearity. 
We give here a brief description of the results of the present paper. 

We first extend the result given in to cover global well-posedness for the problem 11.11 in the 
energy space {Q = H^{M.~^) © H^{W~), coinciding with the form domain of H^) in the subcritical 
case n < 2 and local well-posedness for n > 2. In particular, the equation has two conserved 
quantities, the energy E and the charge (i.e. the L^-norm) M, as in the unperturbed NLS equation. 
They are associated to symmetries of the space of solutions, respectively time translation and phase 
invariance. The free NLS equation has one more conserved quantity, linear momentum, associated to 
space translation; this symmetry is broken by the defect, and correspondingly there is not momentum 
conservation. 

In the main part of the paper we are concerned with the identification of the ground states for 
equation 11.11 and the analysis of their stability in the case of attractive 5' interaction. 
A nonlinear standing wave, or a nonlinear bound state in physical terminology, as in the linear case 
is a solution of the form 

Correspondingly, (j)^^ fulfils the stationary equation 

H^(l)^-\\(l)^\'^''(l)^ = -ujct)^. (1.2) 

We call A the set of stationary states of the equation II. 1[ 
Every member of A has to be a classical (C^(M^)) and square integrable solution of the standard 
NLS to the left and to the right of the singularity. This gives, for the spatial part of the standing 
wave, the only possible forms 

i_ j_ j_ _j_ 

ibA + i)2Mtj2M cosh i^[^y/u:{x — Xl)], x<0 

A 2n (^ + 1) 2Ma;2M cosh t^[iiy/uj{x — X2)]-, X > 0. 

Note that in the introduction we omit, to simplify notation, the dependence of (p from every parameter 
other than the frequency uj. The standing wave solution is represented as a solitary wave of the NLS 
centered at xi on the left of the origin and a solitary wave of the same NLS centered at X2 on the 
right of the origin; the parameters xi, X2 defining the solution are to be chosen in such a way that the 
function (f)^ satisfies the boundary conditions embodied in domain of H^. Eventually, they depend 
on the parameter A, ^, 7, oj which enter the equation. It turns out that there exist two families 
of stationary states, a family Ti whose members respect the symmetry xi = — X2 and a family 
whose members do not enjoy this symmetry. 

It is an important point that the analysis of standing waves by ODE methods has a variational 
counterpart in the fact that the standing waves turn out to be critical points of an action functional. 
The action functional for our problem is defined as 

SM = If,{cP) - ^^iml^^Xl + ^ml = E{<p) + |m(<a), <A e Q, (1.3) 
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where we indicate by F^{(l)) the quadratic form associated to the S' point interaction H^. 
It is easy to see that solutions of (jl.2p are critical point of (jl.3p . The above action is easily seen to 
be unbounded from below for focusing nonlinearities. We show however that the action attains a 
minimum when constrained on the natural constraint (5*^ (</>), i;^) = (the so-called Nehari manifold) 
which obviously contains the set of solutions to Euler-Lagrange equations for S^. In the present 
paper we adhere to the customary mathematical use to call ground states the minimizers of the 
action on the natural constraint. To prove that the constrained minimum problem has a solution we 
exploit a) the boundedness from below of the action on the associated Nehari manifold; b) the fact 
that the Nehari manifold is bounded away from zero; c) classical Brezis-Lieb inequalities showing 
that the limit of a minimizing sequence exists and it is an element of the minimization domain; d) 
finally, minimizers turn out to be elements of D{H^) and satisfy equation (jl.2p . 
For an analysis of the analogous problem in the simpler case of NLS with a 5 interaction see [17j . 
It turns out that ground states do not exist for oj < and that for every u > ^ there is at least 

one ground state. More precisely, for every w G (^, ^%^), there exists a unique (up to a phase 
factor) stationary state. Furthermore, it belongs to the family J^i, so it is symmetric with respect to 
the defect. From such a family, a couple of non symmetric stationary states, i.e., belonging to the 
family J-2, bifurcates in correspondence of the value uj* = 

Before focusing on the issue of determining the stability of the stationary states, we recall that, to 
the ?7(l)-invariance of equation (jl.ip . the appropriate notion of stability for our problem corresponds 
to the so-called orbital stability, which is Lyapunov stability up to symmetries. The stationary state 
(pu) is orbitally stable if given a tubular neighbourhood U (cpui) of the orbit of the ground state (i.e., a 
neighbourhood modulo symmetries), the evolution -(/^(t) is in U{(j)^) for all times if the "0(0) datum 
is near as it needs to the orbit of (j)^^ . Precise definitions are given in section 16.31 

There are two main approaches to determine orbital stability or instability of stationary states: 
the variational or direct method, pioneered in the classic paper |8], and linearization. In our analysis 
we employ the linearization method, studied in a rigorous way first by Weinstein ([32] |33j ) and then 
developed as a general theory for hamiltonian systems by Grillakis, Shatah, and Strauss ([21], |22|). 
"Linearization" in this context is a conventional denomination because the theory includes in fact 
the control of the nonlinear remainder. 

A summary of the essential steps of the method can be sketched as follows. The equation (II. ip 
can be written as a hamiltonian system on a real Hilbert space after separating real and imaginary 
part oi ip = u + iv. Given a stationary state, the linearization of the so constructed hamiltonian 
system at the point cj)^ is described by the second derivative of the action S"{4>u))- Such a quantity is 
identified with a linear operator, whose action is suitably represented with the aid of operators Lj'^ 
and -^2"^ defined through 

S" {(l)u)){u + iv) = Lj'^u + iL^'^^v 

(we refer to 16.21 and subsequent comments for the explicit definition) . Lj'^ and L2 '"^ are easily seen 
to be self-adjoint on their natural domains. 

Then, denote by n{Lj'^) is the number of negative eigenvalues of Lj'^ , and define the function 



1 if £,suM>o 

if ^SM<0 



Now, provided that 

a) the essential spectrum of S'^{(t)ui) is bounded away from zero, 

b) Ker(L2''^) = Span{(/)^}, 
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c) n{Ll'^) = n, 

the stationary state (j)^] is stable if n — p = 0, and unstable if n — p is odd. 

We accurately compute -^Si^{(j)i^) for the ground states of our model, in critical and supercritical 
regime, and prove the occurrence of an exchange of stability between the two subfamilies of J-i. For 
subcritical and critical nonlinearity power /x < 2 and low frequency a; G ^%^), the symmetric 
(odd) states are stable, while crossing the critical frequency w* the symmetric branch of stationary 
states becomes unstable and the two newborn asymmetric states prove stable. After bifurcation 
(w = w* + 0), ground states are stable for not too strong nonlinearities 2 < < //*, where /i* lies 
between 2 and 2.5, and they become unstable for /i > /x*. Finally, in the supercritical regime /i > 2 
and large frequency w ^> a;*, all ground states are unstable. 

We are then in the presence of a pitchfork bifurcation, accompanied by a spontaneous symmetry 
breaking on the stable branches of the bifurcation. The phenomenon of bifurcation of asymmetric 
ground states from a branch of symmetric ones was discovered by Akhmedeev ([3]) in a model of 
propagation of electromagnetic pulses in nonlinear layered media, and then studied in a rigorous 
way in several mathematical papers (see [29] and references therein) with various generalizations 
of the original result, which was concerned with an exactly solvable example. In these works the 
bifurcation is induced by a change of sign in the nonlinearity (i.e., a transition to a defocusing regime) 
in a localized region. Here we have an analogous effect in a different model, in which the bifurcation 
is induced by a strong point defect. 

We remark that no such a phenomenon shows up in the case of a single 5 perturbation of NLS. 
Nevertheless, a situation similar in some respects appears in the case of the NLS with two attractive 
5 interactions separated by a distance, a double 5 - well. This model is studied in [26], where 
an analysis of the model is given by means of dynamical systems techniques, and [T8], where the 
bifurcation is explored in the semiclassical regime. See also [27] for the analogous phenomenon 
with a regular potential of double well type and |30| for the introduction of a related normal form. 
The analogy is a reminiscence of the the fact that the 6' interaction can be approximated in the 
norm resolvent sense by a suitable scaled set of three 5 interactions (see and |141), so it may 
be considered as a singular model of a multiple well. We finally recall that a different definition of 
ground state for a semilinear Schrodinger equation requires that they are minimizers of the total 
energy constrained to have constant mass. The two definitions are related, but not in an obvious 
way, because the constrained action is bounded from below irrespective of the power nonlinearity 
and the constrained energy on the contrary is bounded from below only in the subcritical case n < 2. 
On the other hand, the general GSS theory guarantees that the ground states are local minima of the 
constrained energy if and only if they are stable, and this gives a connection between two definitions. 

The paper is organized as follows: after giving the main notation, we introduce the model by 
recalling few elementary properties of the NLS and of the point interaction (sections 2.1. and 2.2). 
Then we state the variational problem that embodies the search for the ground states and two 
variants of it (section 2.3). In section 3 we establish the well-posedness of the dynamical problem, 
that turns out to be only local in time for strong nonlinearity (/i > 2) and global in time for weak 
nonlinearity (p < 2). In section 4 we prove, by variational techniques, the existence of a ground 
state (theorem 14. ip . while in section 5 we explicitly compute the ground state (proposition 15. ip and 
show the occurrence of a bifurcation with symmetry breaking (theorem 15. 3p . Section 6 is devoted 
to the issue of determining the stability and instability of the ground states: first, in section 6.1, 
we prove the spectral properties a), b), and c) for S'^ evaluated at the ground states (propositions 
16. H 16.31 and l6.4p : then, in section 6.2, we compute ■^Si^{ipi^) and study its sign (proposition 16. 5p : 
finally, in section 6.3 we collect all information and prove stability, instability and the occurrence of 
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a pitchfork bifurcation (proposition 16.91 and theorem I6.1ip . 



1.1 Notation 

For the convenience of the reader, we collect here some notation that will be used through the paper. 
Most symbols will be defined again, when introduced. 

- ujQ := ^ is the frequency of the linear ground state. 

- uj* := ^^^^ is the frequency of bifurcation. 

- The symbol /' denotes the derivative of the function / with respect to the space variable x. Time 
derivative is explicited through the symbol dt- 

- We denote the L^-norm by || • When p = 2 we omit the subscript. The squared H^-norm of / 
is defined as the sum of the squared L^-norms of / and of /'. 

- The following abuse of notation is repeatedly used: 

(f-—e i-oo 
J \^'{x)\Pdx + j \iP'{x)\Pdx, l<p<oo 

- Our framework is the energy space 

Q := H\m+)®H'^{R-) 

ml ■■= \m' + u'f 

- The ordinary scalar product in is denoted by (•, •) and is antilinear on the left factor. We often 
use the duality product in the space Q, denoted by {f,g), where f £ Q* and g £ Q- Again, it is 
antilinear in /. For simplicity, in the bracket (•, •) we often exchange in the bracket the place of the 
element of Q with the place of the element of Q* . 

- The symbol e denotes the sign function. 

2 Preliminaries 

Here we recall some well-known facts. 
2.1 Free stationary NLS 

The set of the solutions to the free (i.e. without point perturbation) stationary Schrodinger equation 
with focusing power nonlinearity 

- ij" - X\ip\'^''ip = -ujip, a;>0, A>0, /i>0 

is given by {4>^,xq € M, w > 0}, where 

(t)^'>{x) = \'^{^JL + l)^uj"^sec\v^[^y/uj{x - Xq)] (2.1) 
Note that they are smooth, exponentially decaying functions. 
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2.2 5' interaction: hamiltonian operator 

The hamiltonian operator describing the so-called 5q interaction is defined on the domain 

D{H.,) = {V' G F'(M\{0}), ^'(0+) = V'(O-) , V(0+) - ^(0-) = -7V''(0+)} (2.2) 
and its action is given by 

H^ij = , X / . 

Note that in the current literature the opposite sign of 7 appears in the definition of H^. This is 
a matter of convention, and we prefer to adopt the present one because we consider the case of 
attractive 6' interaction only. 

An equivalent description of the elements of the operators domain D{H^) is 

V' = (2.3) 

where 

i}± G i/2(M\{0}) ni7i(M) 
are even functions which satisfy the boundary conditions 

^^+(0) - V-(0) = -7^^.(0) = -7V'-(0). 

In this representation, the action of is 

H^Hx) := -x+{x)iIj'1{x) - x-{x)ip'!^{x). 

Let us recall the main spectral properties of H^. The essential spectrum is purely absolutely contin- 
uous: aess{H-y) = [0,-|-oo). 

Concerning the discrete spectrum, if 7 < 0, then ap{H^) = 0; if 7 > 0, then there exists a unique 
eigenvalue, given by (Tp{H-y) = { — ^}. The corresponding normalized eigenfunction is given, V 7 S 
(0,+oo), by 

2\5 _ _2|^| , ^ X 



ip-yix) = — e(x)e 7I I J eix 

Finally, the singular continuous spectrum is empty: (Jsc{H^) = 0. 
In the following we use the notation ujq = 

2.3 (^'-interaction: quadratic form 

It is known (see \l\ that the description of the dynamics generated by can be equivalently 
accomplished by using the quadratic form 

F^(^) := ||^'||2_^-i|V;(0+)-V(0-)P, (2.4) 

where, with an abuse of notation, we denoted ip' := X+V'V + X-V'-i according to the notation 
introduced in Section 1.1. The domain of F^/ is given by (see also Section 2.2 in [IJ) 

D{F^) := H^{W^)®H^{^'). (2.5) 
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We will consider the form domain as the natural energy space for the 6' interaction; it is independent 
of 7 and from now on we will indicate it by Q. If ■0 belongs to the operator domain of a ^'-interaction 
with strength 7, then one has 

:= u'\\'-i\i^'m' . 

which explains the name given to the operator. 

Due to ()2.5p . it is natural to endow the space Q with the norm 

2 l|i ;. ||2 , l|i ;. ||2 ii ;.||2 , |L;/||2 



Q ■- 2^\''y+^^H^R) + 2^\v^-\\miR) - w\\ + WW w ■ 

It is immediately seen that the following proposition holds 

Proposition 2.1. The form is continuous in the topology induced by the norm \\ ■ \\q. 
Proof. In order to estimate the pointwise term in F-y, notice that 

|V(0+)-V(0-)|' < 2|V'(0+)|2 + 2|V'(0-)|2 < 2||V'+||||V;il+2||^-||||V-|| < 211^^111, (2.6) 

where we used the well-known fact that ||</>||^ < 11011211 ^'Ib for any (p £ //^(M). The first term in 
(j2.4p is trivially estimated by the squared Q-norm, so the proof is complete. □ 



Finally, we notice that, for any m > —ojo one has 

F^(^) = {{H^ + m)^,{H^ + m)^ij)-m{'il^,'il^). (2.7) 
This is immediate for ip £ D(H^), then a density argument shows that (|2.7p holds for any ip £ Q. 

2.4 Functionals and variational problems 

We first define the functionals we use through the paper. All of them act on the energy space Q, 
defined in (|2.5p . 

We define the Hamiltonian of the NLS with a point defect as the sum of the linear Hamiltonian of 
the corresponding point interaction and of the nonlinear self-interacting part. 

So, in the particular case of the interaction, or 6' defect, we have for the total energy 

EW = - ^mO+) - ^P{0-r - ^^11^112^2 • (2.8) 

Standard results in the calculus of variations (see for example [10] show that E G C^(Q,M) and the 
Prechet derivative 

E'ii^) = H^^ - X\^|J\^^'^|; £ Q* ^ip £ Q . 
Moreover we define the mass functional (sometimes called charge) 

M{tp) : = 



The mass obviously belongs to C^{Q,M). Both mass and energy are conserved by the flow (see 
Proposition 13. 4p . 

A prominent role in the variational characterization of stationary states is played by the action 
functional 

S^:= E + ^M , (2.9) 
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or explicitly 



i||^'||2 + ^l^||^||2_l 

2"^" 2"^" 27'"' ' " 2/X + 2' 



SM = ^ll^'f + -IIV'IP - -1^(0+) - V(0-)P - TTT^II^II^S' • (2-10) 



Again, S^j G Ci(Q,M) and for every ^ G Q 

SU^P)=H^^P + ui;-X\i;\^^i;GQ* . 

Stationary states V'w satisfy S'^i^Tp^) = . 

Moreover, it is useful to define the so-called Nehari functional 

IM ■. = {S'M)A) =||V'f --|V'(0+)-V(0-)p-A||^g:^ 

(2-11) 

It is immediately seen by its definition that the zero-level set of the Nehari functional contains the 
stationary states associated to the action 5^,. 

The action restricted to the Nehari manifold gives the last auxiliary functional 

:= ^^^^^WntXl = - ■ (2-12) 

In many physical contexts, the search for the ground states can be formulated as follows: 

Problem 2.2. Given m > 0, find the minimum and the minimizers of the functional E in the energy 
space Q under the constraint M = m. 

Nevertheless, in the investigation of orbital stability of stationary states it proves useful to study 
another variational problem, namely 

Problem 2.3. Find the minimum and the non vanishing minimizers of the functional S^^ in the 
energy space Q under the constraint /^^ = 0. 

This is the problem studied in the present paper. It can be rephrased to the issue of characterizing, 
as explicitly as possible, the function 

d(uj) := mf{SM,iP G Q\{0}, 1^ = 0}. (2.13) 

Finally, studying problem 12.31 is equivalent (see Section [4]) to the 

Problem 2.4. Find the minima and the non vanishing minimizers of the functional S in the energy 
space Q under the constraint /(^ < 0. 

Problems 12.21 and 12.31 are related, but not in an obvious way. Of course, when 12.21 has a solution 
ip, via Lagrange multiplier theory it turns out that there exists a real multiplier uj such that E^ip) + 
ujM'(^) = 0, which coincides with S'^{ip) = 0, meaning that ip \s a, stationary point of and by 
definition it belongs to /(^, so it is a solution of l2.3l Nevertheless the two problems are not equivalent 
and a complete analysis will be given elsewhere. See however the remarks at the end of Section 4 
and Section 6. 
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3 Well-posedness and conservation laws 

Here we treat the problem of the existence and uniqueness of the solution to equation (jl.ip . In the 
present paper we are mainly interested in solution lying in the energy space Q and in the operator 
domain D{H^), i.e. weak or strong solutions respectively. Let us stress that it is possible to obtain 
local well-posedness in by proving suitable Strichartz estimates and then following the traditional 
line (see e.g. |10]): this route is followed for a general point interaction and a cubic nonlinearity in 
[1], to which we refer. 

We begin with weak solutions and to this end, instead of equation (jl.ip we study its integral form, 
i.e. ^ 

^(t) = e-'^-'^iPo + iX [ e-'^-'^^-'^\7p{s)\'^f'iP{s)ds, (3.1) 

Jo 

where the initial data ipo belong to Q. 

For an exhaustive treatment of the problem given by a general point interaction at the origin and a 
cubic nonlinearity, see ^j. In particular, we recall that the dual Q* of the energy space Q, i.e. the 
space of the bounded linear functionals on Q, can be represented as 

Q* = //-^(M) eSpan((5(0+),(5(0-)), (3.2) 

where the action of the functionals 5(0ib) on a function ip £ Q reads 

{5{0±),^) = ip{0±). 

As usual, exploiting formula (12. 7p . one can extend the action of to the space Q, with values in 
Q\ by 

{H^ipi,7p2) := {{H^ + m)^iJi,{H^ + m)^ij2) -m{i;i,ij2), (3.3) 

where m > —coq. The continuity of H^tpi as a functional of the space Q is immediately proved by 
Cauchy-Schwarz inequality, (j3.3p and (j2.7p . that together give 

\{H.,i;uiJ2)\ < C||Vi||q||V2||q. (3.4) 
Now we are ready to state the following lemma. 
Lemma 3.1. For any ip £ Q the identity 

^e-'^-'^ij = -iH^e-'^'<^tl: (3.5) 

holds in Q* . 

Proof. The time derivative of the functional e~*^T*^/^ is defined in the weak sense, namely 

\dt I h^oh 

Now, fix ^ in the operator domain D[H^) defined in (|2.2p . Then, 

where we used ()3.3p . 

Then, the result can be extended to ^ £ Q by a density argument, and by (13. 4p we get the 
continuity of the functional -4e~*^^*^ on Q, so the result is proven. □ 



(g-iH^(t+h)^^ .) - (e-^^^V, •) 
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Corollary 3.2. By (j3.5p . one immediately has that the formulation (jl.ip of the Schrddinger equation 
holds in Q* . 

In order to prove a well-posedness result we recall from Section 4 in [1] that the one-dimensional 
Gagliardo-Nirenberg estimates holds for any ip in Q, i.e. 

iiV'iip ^ cu'\\i;hmltK v-gq (3.6) 

where the C > is a positive constant which depends on the index p only. 

Notice that from inequality (|3.6p one immediately obtains the Sobolev-type estimate 



2^+2 < c^IIV'IIq. (3.7) 

Proposition 3.3 (Local well-posedness in Q). Fixed ipo E Q, there exists T > such that equation 
dH]) has a unique solution G C^{[0,T),Q) n C^{[0,T),Q*). 

Moreover, eg. ()3.ip has a maximal solution ip"-'^^^ defined on an interval of the form [0,T*), and the 
following "blow-up alternative" holds: either T* = oo or 

lim ||V''"^^(t)||Q = +00. 

Proof. We denote by the space L°°([0, T), Q), endowed with the norm IIV^IIa' := sup^g^ ||^/^(t)||Q. 
Given ipo G Q, we define the map G : X — > X as 



G<i) ■= e-'^-'-tjjo + I e-^^^(-")|0(s)|2/^</.(s) ds. 
Jo 



Notice first that the nonlinearity preserves the space Q. Indeed, referring to the decomposition ()2.3p . 
since both ip± belong to H^{M.), the functions |V'±PV'± belong to i?^(R) too, and so the energy space 
is preserved. 



Now, estimate (|3.6p with p = oo yields 

ms)\'^Hs)\\Q < cms)\\i^^\ 

so 

\mu <m\Q+c ["^uis^^+Us < Uo\\q+ctu\\i^+\ (3.8) 

Analogously, given (p,^ G Q, 







\\G<p - Geiu ^ CT [uwi^ + mi^) 110 - . (3.9) 

We point out that the constant C appearing in (j3.8p and (|3.9p is independent of 4'i aiid ■C- Now 
let us restrict the map G to elements 4> such that \\(p\\x ^ 2||V'o||q- Prom (13. 8p and (13. 9j) . if T is 
chosen to be strictly less than (SCH^/^oIIq)""^; then G is a contraction of the ball in X of centre zero 
and radius 2||^/;o||Q) and so, by the contraction lemma, there exists a unique solution to ()3.ip in the 
time interval [0, T). By a standard one-step bootstrap argument one has that the solution actually 
belongs to C^{[0,T),Q) and, due to the validity of (jl.lj) in the space Q* (see Corollary 13. 2p . we 
immediately have that the solution ip actually belongs to C^{[0,T),Q) n C^{[0,T),Q*). 

The proof of the existence of a maximal solution is standard, while the blow-up alternative is a 
consequence of the fact that, whenever the Q-norm of the solution is finite, it is possible to extend 
it for a further time by the same contraction argument. □ 
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The next step consists in the proof of the conservation laws. 

Proposition 3.4. For any solution ijj G C^{[0,T),Q) to the equation (|3.1|) . the following conserva- 
tion laws hold at any time t: 



/-oil, E{m) = E{ijo), 
where the symbol E denotes the energy functional introduced in (12. Sp . 

Proof. The conservation of the L^-norm can be immediately obtained using Lemma [3.1l and Corollary 
[321 So, 

(i 2 / , /.N d 



by the self-adjointness of H^. In order to prove the conservation of the energy, notice that {ip{t), H^7p{t)) 
is differentiable as a function of time. Indeed, 

^[{^{t+h),H^^l^it+h))-{m,H,m)] = ( "^^^ "^^^ , H-yi^it + /i)^ + (^H,m, "^^^ 

and then, passing to the limit /i — )• 0, 



m).H^m) = 2Re{-m.H-,4it)) = 21mmt)\''iit),H^m). (3.10) 



where we used the self-adjointness of and corollary 13.21 Furthermore 



l{^{t),\m\'^m) = j^mt)m,nt)m) = {2t,+2)im{\m\'^m,H^m)- (3.11) 



From ()3.10p and (j3.1ip one then obtains 

j,E(m) = l|(.«,).^f,^(,)>-^|»(,),i«*)i^"«*)) = 

and the proposition is proved. □ 
Corollary 3.5. For fJ. < 2, all solutions to ()3.ip are globally defined in time. 

Proof. By estimate (|3.6p with p = oo and conservation of the L^-norm, there exists a constant M, 
that depends on ipQ only, such that 

i?(^o) = E{m) > ^||V'(i)lP-M||^'(t)|| 

Therefore a uniform (in t) bound on HV^'COlP obtained. As a consequence, one has that no blow- 
up in finite time can occur, and therefore, by the blow-up alternative proved in theorem ()3.3p . the 
solution is global in time. □ 

Using a well established analysis which we will not repeat here (see again [1] for a detailed study 
in the case /i = 1)), one can get well posedness for strong solutions of (13. ip . i.e. in the operator 
domain. 
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Theorem 3.6 (Local and global well posedness in the operator domain). For any /i > and initial 
datum TpQ G D(H^) there exists T G (0, +00) such that the equation (jS.ip has a unique solution 
i) e C{[0,T),D{H)) n C^([0,r),L2(M)). Moreover, equation (|^ has a maximal solution ij{t) 
defined on the interval [0,T*). For such a solution, the following alternative holds: either T* = +00 
or 



lim U{t)\\H = +00. 

Moreover, if < fx < 2 and ipQ £ D(H^) the equation (jS.ip has a unique global solution ip G 
C{R,D{H^))nC^{R,L'^{R)). 

Finally, mass and energy are conserved quantities for both weak and strong solutions. 

4 Existence of a ground state 

In this section we show the existence of a solution to problem 12.31 for any u > loq. More precisely, 
we prove the following theorem: 

Theorem 4.1. Let uj > . There exists ^ G Q\{0}, with luii') = 0, that solves problem \2.4\ 
namely 

sui^) = d{oj). 

In order to prove theorem 14.11 we need four preliminary lemmas. In the first lemma we show that 
problems 12.31 and 12.41 share the same solutions. 

Lemma 4.2. For the functionals S^j and S, defined respectively in (j2.9p and (|2.12p . the following 
equalities hold: 

d{uj) := mi{SUi^), G Q\{0}, 1^ = 0} = ■m^{S{^P), V G Q\{0}, Ui') < 0}. (4.1) 

Furthermore, a function cj) G Q\{0} satisfies S{(j)) = d{u)) and Iuj{(p) <0 if and only if Su){(j)) = d{uj) 
and Iui{(j)) = 0. 

Proof. Let </) be a nonzero element of Q such that = 0. Then, by (j2.1ip . Su;{(p) = S{(j)), 

therefore 

M{SM, ^ G Q\{0}, = 0} > inf{5(V), V G Q\{0}, IM < 0}. (4.2) 

On the other hand, let (j) be an element of Q such that Iuji<P) < 0. Defined 

A5i||.*|l2p+»2 

one can directly verify that a{cj}) < 1, I^{a{(t))(j)) = 0, and then, by (|2.1ip . 

S^{a{^)(t>) = ^(a(0)(A) = a(0)2'^+2^(^^ < g^^^^ 

so ^ 

mi{SM. e g\{0}, IM = 0} < inf{S(^)), ^ G Q\{0}, 1^ < 0}- (4.4) 

From (ji?2]) and identity ([i?T]) is proven. 
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Prom ()4.ip . it is obvious that, if (p minimizes S^j on the set /(^ = 0, then it minimizes S on the 
set /oj < too. ^ ^ 

Suppose now that £'((/>) = d{uj) (then minimizes on the set < 0) and Iuji4') < 0. Defining 
a{(f)) like in one finds a{<f>) < 1, Iui{a{<f>)(f)) = 0, and 5a;(a((^)</)) = a{(f))'^^'+'^S{<f>) < d{oj) again, 
that contradicts the definition of d{uj). So the lemma is proven. □ 

Lemma 4.3. The quantity d{uj) defined in ()2.13p is strictly positive. 

Proof. First notice that, for any a > 0, 

|V^(0+)-V(0-)P < 2|V'(0+)|2 + 2|V(0-)|2 < 2||^+||||V'^||+2||^ 



where we used decomposition (j2.3p . estimate (13. 6p with p = oo, and Cauchy-Schwarz inequahty. 
Therefore, 



F,{,l:)+u:Ur > [l-^^j\Wf+[u-^-^^Uf. (4.6) 



Since uj > we can fix the parameter a in such a way that 

2 -fuj 
- < a < —, 
7 2 

so it is proven that 

F^(V^) +^11^^11' > Cll^ll^. (4.7) 

Then, by the estimate ([37 



IM > C7||V||^-A||Vll2;:+2 > C^lllV'lli+2--ll<.ll2/.+2- 



|2/i+2 ^ r*, IU/,I|2 _ _ \IU/,l|2/'+2 
12(1+2 

It appears that, if luji^P) < 0, then either ^ = or 



>[-) > 0. 

Therefore, since in problem 12.41 we are looking for a non vanishing minimizer, and owing to the fact 
that on the Nehari manifold S^^ = S, it must be d{u}) > 0. 

□ 

In the third lemma we consider a pair of functionals 5^, I^, that correspond to the functionals 
Suj, lui in the absence of the point interaction: 

Sim = ^iiv^'f + 111^^11^-^11^11^X2 

I'M = U'f + u;Uf-xMl>;Xl. 

Lemma 4.4. The set of the minimizers of the functional among the functions in Q\{0}, satisfying 
the constraint = 0, is given by 

{e^'x+<Pl e^\-0°, 0G[O,27r)}, 

where the function cp^ has been defined in (j2.ip . 
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Proof. First notice that, reasoning like in the proof of lemma [^^21 one can prove that the search for 
the minimizers of among the nonzero functions in Q that satisfy = 0, is equivalent to the 
search for the minimizers of 5 among the nonzero functions in Q that satisfy < 0. 
Let us define the real function of a real variable 

d°(a;) := inf{5(V), ^P G Q\{0}, < 0}. 

Proceeding like in lemma HTSl one can show that dP{(jj) > 0. 

Besides, we recall that (j)^ minimizes the functional S among all functions in i/^(M)\{0} such 
that = 0. Now, making resort to the representation (j2.3p . let us consider a generic function of 
(5\{0} supported on M+, call it with G H^(M.) and even, and suppose that I^{x+ip+) < 0. 

One immediately has 

so x+</'[!) is a minimizer of S among the functions of (5\{0}, supported on and satisfying < 0. 
Notice that the equal sign holds if and only if -0+ = 02,. Otherwise, t/j^ would not belong to the 
family (|2.1|) . nevertheless, as S^{^p+) = S^{(l)^), it would be a minimizer of among the nonzero 
functions in H^(M.) that satisfy = 0, which is impossible. 

Thus, for any function ip G Q\{0}, with ijj = x+'i'+^X-'^-j ^ 0) the following alternative 

holds: either < and so 

Siij) > > SU+<Pt), (4.8) 

or ll^{x-ip-) < and so 

S{i,) > su-i^-) > SU-4>'J, (4.9) 

and the equality in the last step of (j4.8p and (|4.9p holds if and only if \ip+\ = (p^, or = cp^, 
respectively. 

Taking into account the [/(!)— symmetry of the problem, the proof is complete. □ 
Lemma 4.5. For the infimum of problem \2.4\ the following inequality holds 

dioj) < S{e'\+<lt) = ^(^^y J+-^ j\l-u')Uu. (4.10) 
Proof. We notice that the last identity in (j4.10p can be obtained by direct computation. Furthermore, 

Ux+h) = Ux+<Po)--y — ^ — j ^ ~i\ — A — ) ^ °- 

Following the proof of lemma 14.21 we define 

F,{x+ct>l)+o,\\x+(ttf]^- _ I, 1 2|<A°(0+)|2\^ _ / ^^ \ 



a : 



>^^^\\xmI'!::2 V ^^mlXlJ \ 7(M+i)^^/;(i-n2)^rf^. 

Therefore, /^('^/t+'/'o ) = aiid 



'^u}\\2fj.+2 

< 1. 
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and, since a < 1, we get 



1 



and the proof is complete. □ 
Now we can prove theorem 14. 1[ 

Proof. Let {ipn} be a minimizing sequence for the functional S on the set I^j < 0. We show that 
there exists a subsequence of {^pn} that converges weakly in Q. First, notice that HV'nllQ is bounded. 
Indeed, the sequence ||V'n||2^t+2 is bounded as it converges. Furthermore, by the lower boundedness 
of the form F^, and recalling that Iu){4'n) < 0, we have 

< < F^(V'n)+a;||^„f < XUnWZXl < C, 

so llV'nII < C. 

Then, using Iuj{ipn) < again, by the decomposition (12. 3p and estimate (13. 6p . 



^if < X\\^l;X'utl-^Unf + -\MO+)-MO-)\' < C + -{\^Pn,+m' + \i^n,-m') 



7 7 



< C + C||V^„,+ |||K ,||+C||Vn,-|||K,-|| < C + C{-Unf + eWJ^ 



(4.11) 

Choosing e sufficiently small, we obtain that HV'nlP bounded, so the sequence {ipn} is bounded in 
Q, and then, by Banach-Alaoglu theorem, there exists a converging subsequence, that we call {V'n} 
again, in the weak topology of Q. We call ipoo the weak limit of the sequence {ipn}- 

We prove that ipoo 7^ 0. To this aim, we show, first, that the sequences {V'n(0±)} converge to 
V'oo(0±), and, second, that lim„_^oo lioi'^n) = 0. 

Let us define the functions '■p±{x) := x±{x)^'^^ ■ Then, integrating by parts, 

+00 f+ca 



{ip+,1pn)Q = e ''Vn(x)dx- / 6 ""ip'^ix) dx = Vn(0+). 

Jo Jo 
Analogously, {ip^,ipn)Q = ^n(0— ). Therefore, by weak convergence, 

^pn{0±) = {^±,^Pn)Q {^±,'>Poc)q = M0±), (4.12) 

and the first preliminary claim is proven. We prove the second claim by contradiction, i.e., supposing 
that Iuj{ipn) — is false. Then, there must be a subsequence of {ipn}, denoted by {ipn} too, such 
that 

lim I^{iPn) = -/3 < 0. 

n— >oo 

We define the sequence Cn '■= ^ni^n-, with 



T^n — ^ ^^1 < i. 
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Since 



we obtain 



lim Vn = lini 



1 + 



^IIV'n|l2^^+2, 



1 

2/^ 



2{fi + l)d(w) 



1 

2n 



< 1, 



lim 5(C„) = Z^'^+'5(Vn) < SiiJn) 



and, since Iu}{Cn) = 0, it follows that the assumption that {V'n} is a minimizing sequence is false. 
Therefore, it must be 

lim /^(V^„) = 0. (4.13) 

To prove that ipoo 7^ we proceed by contradiction again. Assume that ij^ao = 0. Define the 
sequence 7/„ := Pn'4'n with 



Pn := 



[\Wnr+oo\ 



2(j 



(4.14) 



Using and (liT^ we obtain 



lim p„ 

n— >-oo 



lim 

n— >oo 



1 + 



/^(V'n)+7~'IV'n(0+)-V'n(0-)|^ 
■^IIV'n|l2jl+2 



2n 



and therefore 



0, 



(4.15) 



lim S{r,n) = lim p^/^+^^^^j ^ ^(^), 

71— ^-OO 71— >CXD 

Moreover, owing to definition (I4.14p . 

I'M = I'Mn) = +a;||V'nf -Ap^^llVnlll:^ 

SO, due to Lemma 14.41 

d{u) > 5°(x+</)o). 
On the other hand, by Lemma 14.51 we conclude 

d{uj) < Six+M < S{r]n). 

that contradicts (I4.15p . So the hypothesis V'oo = cannot hold. 

Now we prove that /aj(V'oo) < 0. To this purpose, we follow the last lines in the proof of proposition 
2 in |17j . First, we recall an inequality due to Brezis and Lieb (jZ]): if tin converges to Uoo weakly in 
LP, then 



IPnIlp - iFn - -Uoollp - Ihoollp > 0, Vl<p<00. 

First, we notice that if n„ = ipn and p = 2p + 2, then (I4.16P yields 

S{tpn) - S{lpn - Ipoo) - S{Tpoo) > 0. 



(4.16) 



(4.17) 



Further applying (I4.16P to the sequence {V'n} and to the sequence {ip'n} with p = 2, and using (j4.12p . 
yields 



0. 



(4.18) 
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Suppose /^(V'oo) > 0. Then, by (ITOjl and (ITOD . 

lim I^ilpn - Ipoc) = lira I^{Tpn) - Iu;{tPoo) = -luiilpoo) < 0. 
n— >oo n— >oo 

Choose n such that luji^Pn — V'oo) < for any n > n. Then, by definition of d(u}) we have 

d{uj) < S{iPn-^oo), yn>n, (4.19) 
but, on the other liand, ipoo 7^ implies S{ipoo) > 0, and, together with (j4.17p . this yields 

lim S{ipn-ipoo) = lim Siipn) - S{ipoo) = - S'(V'oo) < d{u;), 

that contradicts (|4.19|) . and so it must be Iui{ipoo) < 0. As a consequence, by definition of d{uj), 

5'(^oo) > d{uj). 

Now, since V'oo is the weak limit of {ipn} in L?^^"^-, we must have 
which implies 

S(^oo) = d{u:), (4.20) 

and so ipoo is a solution to the minimization problem [231 and therefore, to the minimization problem 
12.31 too. The proof is complete. 

□ 

Corollary 4.6 (Strong convergence). If a minimizing sequence {ipn} for the problem \2.4\ converges 
weakly in Q, then it converges strongly in Q. 

Proof. Formulas (j4.16p and (j4.20p prove that {ipn} converges strongly in L^'^"^^. As a consequence, 



A* + 1 /U + 1 



and by ()4.7p this is equivalent to the strong convergence in Q. □ 

We end this section by adding some remarks on the variational problem 12.21 i.e. to minimize the 

energy at fixed norm. 

To be precise, let us define the manifold 

Tm = {ip(^Q:Uf =m} (4.21) 

and 

-£n. = \ni{E{i,) I V G r^} • (4.22) 

In the following results it is shown that in the supercritical regime the constrained energy is un- 
bounded from below and in the subcritical regime its infimum is finite and negative. Moreover the 
energy is controlled from below by || • ||q norm. 
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Lemma 4.7. (Behaviour of the constrained energy). 

1) Let fi > 2; then S^n = +c« and the energy E is unbounded from below in r^; 

2) Let fj, < 2; then < < +oo; moreover there exist positive and finite constants Ci, C2 
(depending on ij,,j,m) such that 

E{ij)>CMl-C2 VVGr„; (4.23) 

3) Let /i = 2; then, there exists m* > such that for m < m* inequality (|4.23p holds true. 
Proof. To show 1) let us consider the trial function 

/m 

— e 



^(a,x) = 

(27rcj2) 

A direct calculation shows that $(ct, x) € , and 

E{^a,x)) = ^\\^'{l,-)f- 



1^(1, 



I2/X+2 



2;u + 2' 

This proves that for fj, > 2 the energy is unbounded from below. Moreover, for fi < 2 and a big 
enough, E{^{a, •)) < 0. 

Now, let us show the bound ()4.23p . The nonlinear term and the point interaction term in the energy 
are dominated by the kinetic energy. Let us consider first the nonlinear term. 
Gagliardo-Nirenberg estimate ()3.6p jointly with the condition ip € Tm give 



2^+2 



With the use of the classical elementary inequality 



2\1 



X'' y 

xy < 1 

P Q 



one obtains, for any e > 
* = Cm'' 



1 1 

- + - = 1 

V q 



< Cm'' 



/,/||2, 



2 



+ 



e 



2 

2-/1 



(4.24) 



from which it follows 



2 

2^l+2 



r,l2<CieU'\\Q + C2 



(4.25) 

In an analogous way one can treat the point interaction part of the energy. Again taking in account 
that ijj G Tm and by use of Sobolev embedding in one dimension and elementary inequalities, one 
has 

1 1 r ii„/,ii2i 1 

12 +S 



i|V(o+)-V'(o-)p<i 

7 7 



1 

< -e 

7 



(4.26) 



Collecting the estimates for the nonlinear part and for the point interaction part of the energy and 
choosing e small enough one gets points 2). 

To prove 3), it suffices to notice that from (13. 6p one immediately has 



;,'i|2 



for any ip S L^, and the proof is complete. 



□ 
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Note that the constrained action attains its minimum for every positive value of fi, at variance 
with the constrained energy, which is unbounded from below for /i > 2. Even for < 2 it is not 
guaranteed that the energy constrained on has a minimum, i.e. that there exist a solution to the 
variational problem 

-£m = m.m{E{Tp) I V G ^m} ■ 

An analysis of this problem for NLS with point interactions will be given in [2]. However, let us note 
that if a minimum exists at ipm, and the constraint is regular at Vmj there exists a Lagrange 
multiplier Am such that E'{7pm) + A^M' (ipm) = 0. This means that V'm is a stationary point for S^^ 
with CO = Ajn and so iIj^. S Ii^. 



5 Identification of the ground state: bifurcation 

Proposition 5.1. Any solution ip to Problem \2.3\ has the form 

where the functions (j)^ have been defined in (12. ip . 6 can he arhitarily chosen, and the couple (xi,X2) 
is determined in the following way: denoted ti = tanh{fl^/uJ\xi\) , then (ti,t2) solves the system 

U U Z ■ (5-2) 



Proof. Consider the functional Joj^i/ = + i^I^, with i' a Lagrange multiplier. Any solution of 
problem 12.31 must then be a stationary point for Juj^u- Let ip be one of such solutions. Then, 

= = 0. 

Furthermore, 

and therefore J^^^^i'^)'^ — ~'^^f-M\'^\\2^+2- Thus, for nontrivial solutions it must be = 0. We 
conclude that any non vanishing minimizer ijj of the functional J^^,y must fulfil S'{'ip) = 0, i.e. 

S'Mri = 0, Vr? G Q. (5.3) 

Applying ()5.3p first to r], then to = —ir], and summing the two expressions, we find 

B^{^,7])-X{\iP\^f'iP,7])+uj{^P,rj) = 0, (5.4) 

where we used the shorthand notation 

1 



i?^(V,r?) := (V',r?')--(V'(0+)-V'(0-))(r?(0+)-77(0-)). 

7 

So, from ()5.4p the following estimate holds. 

\B^{i^,v)\ < AIIVII^IIV'IIII^II < C^^ll^ll, Vr/eQ. (5.5) 
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Notice that, letting rj vary among the functions vanishing in a neighbourhood of zero, we conclude 
from (jS.Sp that ip G /7^(M\{0}). Thus, for a generic r] £ Q a straightforward computation gives 



^^(0+) - '0(0- 



B^{i^,v) = {-^",v)-{v{0+)-v{0-)) y' ' ' +^'{0-)j, (5.6) 

where we used the notation tp" := x+V'++X-V'" • So, from (jS.Sp and (|5.6p . we conclude that tp belongs 
to the domain D{H^) (see definition (j2.2p ). As a consequence, the function H^ip — Aj-i/'P'^V' + ^'4' 
belongs to ^^(M). 
Furthermore, from (|5T 



and, since Q is dense in L^(M), 

H^^p - X\'p\'^^'^ + uj^ = inL2(R). (5.7) 

Since ip lies in the domain of H^, equation (j5.7p can be rewritten as 

-V'"-A|'0|2a'V' + w'0 = 0, V'GF2(M\{0}) 

^-'(0+) = ^'(0-) . (5.8) 

V'(0+)-V(0-) = -7V''(0+) 

Consider first the case of a real V'- Then, the first equation can be interpreted as the law of motion of 
a point particle with unitary mass, moving on the line under the action of the double-well potential 
^(^) ~ 2/i+2 ••''^^"^^ ~ standard methods of classical mechanics (see e.g. [19]) one immediately 

sees that the only solutions that vanish at infinity correspond to the zero-energy orbits, whose shape is 
given by (j2.ip . where xq is a free parameter that, in the mechanical problem, embodies the invariance 
under time translation. 

Consider now the possibility of complex solutions. Writing ip{x) = e^^^^'^p{x), the first equation in 
dSSD yields 



-P 



-2ie'p' - xp^^'^^ + {uj + e")p = 0, 



thus, in order to make the imaginary part vanish, either p' or 6' must be identically equal to zero. If 
p' = 0, then ^p either vanishes or is not an element of L'^(R). So it must be 6' = 0, and since M\{0} 
is not connected, one can choose a value for the phase in the positive halfline and another value in 
the negative halfline. One then obtains that all possible solutions to (15. 8p must be given by 



= \ , (5.9) 



e'^^^cpZ'ix), x<0 
gie2 0^2(^)^ x>0 

where xi and X2 are to be chosen in order to satisfy the matching conditions at zero. 
We remark that, among the functions in (15. 9p . once flxed xi and X2 the minimum of 5"^^ is accom- 
plished if the condition e*^^ = — e*^^ is fulfilled. Indeed, it is clear that such a condition minimizes 
the quantity — 27~^|^(0-|-) — ip{0—)\'^, while the other terms in the definition (12. 9p of the functional 
Soj are the same. This explains the phase factor in (|5.ip . 

Owing to the phase invariance of the problem, without losing generality we can choose 6i = vr, 
02 = 0, so the matching conditions in ()2.2p yield the following system for the unknowns xi, X2, and 
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oj: 

ta,nh{fj,y^xi) . tanh(/tY^a;2) r> 

1 1 1 — U 

coshi^ (fiy/Uxi) coshM (^y^x2) 

(5.10) 

T 1 T — 1 • 

cosh (/i\/iU^2) cosh'' (/^-y/aJxi) cosh^* (/xy'aJxi) 

By the first equation of system (|5.10p . xi and X2 must have opposite sign. Furthermore, the second 
equation gives xi > 0. So it is proven that X2 < < xi. 

Denoting tj = tanh(/i-y/a;|xj|), and exploiting elementary relations between hyperbolic functions, 
system (j5.10p gives (j5.2p and the proof is complete. 

□ 

Before explicitly showing the solutions to the problem ()5.ip , ()5.2p , we prove a preliminary lemma. 



Lemma 5.2. For any n > 0, a > 2 J there exists a unique x G l] such that 



(ax-l)2M+2 +^ ^ 
Proof. Let us denote t(;(x) = '^^^^lyi^J+^i^'^ + x^ — 1. First, notice that w (|) = 0. Furthermore, 

w'(-] = ^(4(/x + l)-/xa2) 
\ a / a 



as a > 2^ -pp^. Therefore ?i;(x) < in some right neighbourhood of ^. On the other hand, 

w{i) > 0, so the set H whose element are the zeroes of w in (f,l], is not empty. Let us denote 
xi := minH. Then, since w is regular in (^,+oo), it must be either w{x) < or w{x) > in 
some right neighbourhood of xi. In the first case, xi is a local maximum for w. Besides, since 
w{l) > 0, there exists X2 > xi such that w{x2) = 0. As a consequence, there are two local minima 
^ (a'^i)' 2/2 £ {xitX2)- Owing to the mean value lemma, there exist three points zi, Z2 and 
23, lying respectively in a neighbourhood of yi, xi and X2, such that w"{zi) > 0, w"{z2) < and 
w"{z3) > 0. Owing to the mean value lemma again, there exist si £ (21,22) and S2 G (22,23) such 
that u;"'(si) < and w"'{s2) > 0. From the explicit expression 

w"'{x) = [-4/ua3(a2 - 1)(2^ + 1)(^ + 1)2;^ + 4a2(/i + 1)(2/U + l){2fia'^ + 3)x 

+ 40(5^^0^ - 2fi^a^ + fi + 3fia^ - l)]{ax - 1)-2m-5 ^^'^^^ 

it is clear that w"'{x) < for large x.It follows that w'" undergoes at least two changes of sign in 
the interval (^,+c«), but the expression (15. lip shows that in the interval (i,+oo) there is a single 
change of sign only. As a consequence, our starting assumption is false and it must be w{x) > in 
some neighbourhood of xi. 

Let us suppose that there is a point X2 > xi such that w{x2) = 0. Following the same reasoning as 
before, we conclude that w'" must change sign at least twice in (^,+00), that contradicts (15. lip . 
As a consequence, there is only one zero (i.e. xi) of w in (^, +00), so the lemma is proven. □ 

Theorem 5.3. IJooq < oj < u>* , then the solutions to Problem \2.3\ are given by ip^' ^'^ (see definition 
fO)). with e £R and 

,= -i^log2^, (5,12) 
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If 00 > uj* , then the solutions to Problem \2.3\ are given by ^H}' and tp^' , with G M and 



1 + 



y2 = - — ^log 



l+t2 



l-t2 



(5.13) 



1-ti 

where the couple (ti,t2), with ti < t2, solves the system (|5.2p . 
Proof. The function 

/(t) := ^2/^ - t2/^+2 (5.14) 

vanishes at the points and 1, and is strictly positive in the interval (0,1). Furthermore, in the 
interval (0, 1) its only stationary point is i := ^J'-^j, where the function / has a local maximum and 

takes the value m := ■ 

As a consequence, given a > 0, the system 

a = f{ti) = f{t2) (5.15) 

in the unknowns ti and t2, has no solutions for a > m, the unique solution ti = t2 = t for a = m, 
and, imposing ti < t2, three solutions for < a < m: indeed, there exists a unique couple ti,t2, 
with ti £ [0,t), t2 G (i, 1], such that /(ti) = 7(^2) = cl- Therefore, the three couples (ti,ti) (^2,^2)1 
(ti,t2), solve (I^T^ . 

So the set of the solutions to the first equation in (j5.2p with ti < ^2 consists of the union of 



Ti := {0<ti = t2< 1} 

and 

T2 := {(ti,t2), < ti < t < t2 < 1, /(ti) = f{t2)}. (5.16) 

Due to the regularity of /, r2 is a regular curve (see Figure 1). 

We consider the second equation in (15. 2p . Varying the parameter cj, it describes a family of 
hyperboles in the plane (ti,t2); whose intersections with Ti and T2 provide the required solutions to 
the system (j5.2p . 

First, observe that 

rnin (tri+t^-i) = 2, (5.17) 

tl,t2G-'l 

and such a minimum is attained at ti = t2 = 1- 
Second, we claim that 



inf (tr' + i2') =2a/^^, (5.18) 

and such a value is attained at ti = ^2 = t. To show this, we use the Lagrange multiplier method, 
and find that any stationary point of the function t^^ + tg ^ constrained on f{ti) = f{t2) must satisfy 
the condition 

tlfih) = -4f'{t2), ti and t2 + t. 

Let us define g{t) := t^f'(€). Notice that (7 > in (0, f), and (7 < in (t, 1]. Therefore, the condition 
5(^1) = —g{t'i) with < ti < t < t2 < 1 is equivalent to g^{t\) = 5^(*2), < ^2- 
Observe that 



Hi) = / #5'(s)^is = / {^s\r[s)f^2s'f\s)f"{s))ds. (5.19) 
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Figure 1: The system (j5.2p . Full lines represent the solutions to the first equation, while dashed 
lines represent the family of hyperboles given by the second equation. The concavity of the curve 
represented T2 refers to a sufficiently large value of fi (for instance, fi > 1/2). For fi near zero, T2 
can exhibit some changes in convexity, but the result on bifurcation and on the number of solutions 
still holds. 



In the interval [0, t\ the function / is monotone, so it is possible to perform the change of variable 
t^y = fit). Then 

rfih) 

g\ti) = / v{{r{y))dy (5.20) 



where we introduced the function 

p{t) := Aff'{t) + 2t^f"{t) = (8;u2 + 4/x)t2/^+2 _ (g^2 ^ 20// + 12)^2^+4 (5.21) 

and the function r, that is the inverse of / in the interval [0,t|. 

Analogously, exploiting the monotonicity of / in the interval [t, 1] , one can change the variable 
in the integral in (j5.19p and obtain 

rfit2) 

9\t2) = / p{a{y))dy (5.22) 



where the function a is the inverse of / in the interval [i, 1]. 



The function p is non negative in 0, y 2^^+5/1+3 ' vanishes at the endpoints of the same 

interval. Its only stationary point in the interval (0, 1) is located at t := 2^:^+7^1 +6 ^'^^ ^ local 
maximum, where 

In the interval (t, 1] the function p is negative and monotonically decreasing. As a consequence, from 
([021) one gets 

g\t2) = / \p{a{y))\dy. (5.24) 

Jf(t2) 
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Besides, since 



Pit) 



8 



+ l)M+i • 



(5.25) 



one immediately has t > J 2 -tj^^^^ > t- 




Comparing equations ()5.23p and ()5.25p . p (t) < \p{t)\ if and only if 




which holds for any value of fj,. Therefore, max^gjo,*] = 
As a consequence, if si < t < S2, then 



\P{S1)\ < \p{t)\ < \P{S2)\ 



and recalling that ti and ^2 are defined to fulfil /(ti) 



f{t2), 



g\ti) = 2 / p{T{y))dy < 2 



(5.26) 



where in the last identity we used ()5.24p . 

It follows that the only point where g^{ti) = ^^(^2) is given by ^1=^2 = so there are no stationary 
points of the function t^"^ + on the set T2. By comparison with the endpoints (0,1) and (1,0) 
one immediately has that ti = t2 = t corresponds to a minimum, so by direct computation our claim 
(j5.18p is proved. 

As a consequence, from (15.17P and (I5.18P we have that: 

• if oj < ujQ, then the system (15. 2p has no solutions; 

• if ojQ < uj < oj* , then the only solution to (]5.2p lies in Ti and reads = t2 = 2^ ' 

• if a; > w*, then the system (j5.2p exhibits three solutions: the first one lies in Ti and is given 



Furthermore, in the set T2 consider the region ti < < t2', expressing ti from the second 
equation of (j5.2p and plugging it into the first one, one obtains the equation w{t2) = 0, with 
w defined as in the proof of lemma 15.21 and a = 7v^- By virtue of the same lemma, such 
equation has a unique solution ^ in the considered interval. The unique value of ti such that 
(ti,t2) is a solution to (15. 2p is given by (i*)"^ = 7\/w — (^2)^^- 

Due to the symmetry of (j5.2p under exchange of ti and t2, the third and last solution is given 



Any solution to ()5.2p singles out a stationary point of the functional S^j on the Nehari manifold, 
that is unique up to multiplication by a phase. Obviously, the value of on such functions is 
independent of the phase. 

Defining y like in (j5.12p . yi and y2 like in (j5.13p . and owing to (j5.ip . we conclude that: 

• if ujQ < u) < uj* , then the only stationary point (up to a phase) for the functional S^j is given by 
^y,-yfi ^ Due to its uniqueness, it must be the minimizer for S^j whose existence is established 
by theorem 14.11 The explicit expression for y given in ()5.12p is found imposing ti = t2 in the 
second equation of (|5.2p . 



by h = t2 = 



by {t*,q). 
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• For > w* two further solutions appear. Keeping into account the signs of xi and X2 estab- 
hshed in proposition ()5.ip . the two related families of solutions can be denoted by 
^y2,-s/i,6'^ with ui and y2 positive numbers. Obviously, the functional takes the same value 
on them. In order to establish which stationary point is the minimizer we must compare 

su^l'^-'^"') with sui^a^-^^')- 

Let us proceed with such a comparison. From (j2.12|) we know that the functional S^j reduces to S 
when evaluated on stationary states. We have 



and 



Introducing the function ip{t) : 



2A^ 



-+- - 

UJ'^ ^ (/U + 1) 

2Am 

_ t 



[l - dt 



2 



2 

7^ 



dt 



t2 



we obtain 



tj2^''(/i+ 1)" 



We define the function 
thus 



2Am 

CJ^^M (/i + 1) 



[l-t^)T^dt 



{l-f)^dt 



2 

(1 

2 " 



7^^ 



2Am 



{l-f)^dt 



-1 



t2 



[l - t^)^ dt 



q{t) :-- 



fit) 



(1 - Z^2)m du, 



W2 + M (/i + 1)/ 

2A^ 

+ + 1)/ 

2A^ 



il-t^)^ dt 



[I - t^)~i^ dt 



Since 



= {i-t^y-{i-^^{t)y 

the stationary points of q must solve the equation 



t\l-t^)-.-{l-^\t))T^^\t) 



t^ ' 



0, 



(5.27) 



(5.28) 



which is equivalent to the first equation of (|5.2p in the unknowns (t,— 99(4)). Furthermore, from 

(I527D, 

-<f~^{t) = 7V^. (5.29) 



Thus, from ()5.28p and ()5.29p we conclude that the couple {t,—ip{t)) solves the system (15. 2p . This 
implies that the functional S, restricted to functions of the kind (j5.ip . has stationary points only in 
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correspondence with the solutions of (j5.2p . In other words, for > uj* there are three stationary 
points for q, and they coincide with y, 2/1,2/2 • 
A straightforward computation provides 



which is positive if 00 > uj*. Then, is a minimum for q, and since q is regular and there are no 
other stationary points, it must be 

Thus we conclude that V'oj^'"^^'^ and -0^^'"^^'^ are the minimizers, and the proof is complete. □ 

We end this section showing that the branch of nonlinear standing waves bifurcates from the 
trivial (vanishing) stationary state for u) > ojq = ^ . 

Proposition 5.4. Let us consider the branch of symmetric standing waves 
The following relations hold true 

1- II iv.—v.9\\ T II iv.—v.8\, 

\H2{m\{o}) 



lim IIC'-^'^II = hm ||C-^''lb= lim Mr^'^W H2<m{n^^ = 



UJ—^UJQ UJ—^UJQ UJ—^UJQ 

Proof. The result immediately follows by observing that the solution (^~^) tv^) ^° dS^Z]) tends to 
(1, 1) as uj ^ ujQ + 0. So, in the same limit y tends to +00, and the result immediately follows from 
the explicit expression ( (|5.9p . ()5.12p ) for ■0^'"^'"- □ 

6 Stability and instability of the ground states 

6.1 The linearized evolution around a stationary solution 

We study the second variation of the functional S^i- Indeed, for any w £ Q 

-^Su;{ip + sw) = Re{H^^P + ujiP-X\iJj\'^f'iP,w), 

which is often referred to as the fact that 5'^[V'] = H-ytp + w^/; — AjV'P'^V' i^i a weak sense. Therefore, 
linearizing the equation (jl.ip around the real function can be made by computing the second 
variation of the functional S^. Given w £ Q, with u = Retf and v = Imw one has 

1 r 

—^SUi^ + sw) = \\w'f + u\\wf - -\w{0+) -w{0-)\'^ - X / I w{x)\'^\il^{x)\'^^dx 
ds^ 7 Jm. 



2fiX / u^{x)\'ip{x)\'^''dx 



lln'f +a;||nf --|u(0+)-n(0-)|2-A(2// + l) f u{xf\i;{x)\'^f'dx 
7 Jr 

+ \\v'f + uj\\vf --\v{0+) -v{0-)\'^ - X [ v{xf\i;{x)\'^^'dx. 
7 Jr 



(6.1) 
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Defining the operators 



L^v =H^v + uJV-\\il)\^^'v 
on the domain D{Ll''^) = D{Ll''^) = D{H^) (see we get 



(6.2) 



^S^iiP + sw) = (n, Ll^^u) + [v, Ll'^'v). (6.3) 

Now we derive the general spectral properties of the operators L\'^ and Lg''^, needed to prove 
stability or instability of the stationary states. 

It is easy to show that L'\'^ and ^2"^ are self-adjoint operators in L^(M). In fact they are abstract 
Schrodinger operators of the form {H^ + w) + Vi{x) , where the perturbation Vi{x) = Cj|V'P^(x) is 
given by a bounded and rapidly decaying function, and ci = 2^ + 1 , C2 = 1 . Let us consider the 
couple of operators ^^2"^ and —-^ + uj — \\ip\'^^^. Both are self-adjoint extensions of the same closed 
symmetric operator with defect indices (2,2); so their resolvents differ for a finite rank operator. 
As a consequence, thanks to the Weyl's theorem (see [31], Theorem XIII. 4), the essential spectra of 
L]''^ and + oj - AjV'P^ coincide. Moreover, (Te{-^ + uj - AjV'P'') = o-e(-^ + a;) = [w + 00), 
because the potential V2 is (— ^ + a;)-compact. The same reasoning holds for the operator L^'*^, so 
we can conclude 

aeiLj''^) = ae{Ll''^) = [a;, +00). 

Moreover, the fact that Lj'^ and -^2"^ are symmetric and relatively compact perturbations of the 
self-adjoint nonnegative operator + oj allows to conclude (see for example [24J , Theorem 6.32) 
that the possible discrete spectrum is finite or accumulates at the border of the essential spectrum, 
which in our case is positive. So the negative spectrum is finite. 
We will often use the previous remarks without repeating the argument. 

We need more detailed spectral information on the operators Lj'^ and Lg'"^, in particular concerning 
the number of negative eigenvalues. A standard technique to deal with this sort of problems in the 
case of operators with domain elements which are regular enough (typically Schrodinger operator 
with a smooth enough potential) relies on the Sturm oscillation theorem which relates the number 
of nodes of an eigenfunction to the ordering of the corresponding eigenvalue. So, if ip is positive, 
then it coincides with the first eigenfunction, which is simple and corresponds to the ground state. 
This reasoning is not applicable in our case, due to the singular character of H^, with possibly 
discontinuous domain elements. By the way, the problem is not completely settled neither for the 
case of the milder 6 interaction, so we give an independent proof of the relevant spectral properties 
for this case too. 

The results are based on a generalization of a ground state transformation for the singular operator 

H-y. 

Proposition 6.1. Let e'^^^ip{x) be a stationary solution to problem (j3.ip with ip(0-\-)ip(0—) < 0. 
Then, for the operator L'l^'^ defined in (16. 2p . the following statements hold: 
a) Ker L^'^ = Spanj'i/'} 

> 

Proof. Along this proof we denote the operator Ll^'^ by L2 . Proceeding like in the proof of Proposition 
15. H up to formula (15.70 . one obtains that the function -i/; is regular in M\{0} and fulfils the boundary 
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conditions defined by the (^'-interaction. So, by the equation for the stationary states (j5.7p again, 
one immediately verifies that L2V' = and point a) is proven. 

To prove b), notice that for any (f) £ D{L2) the following identity holds at any point x ^ 



1 d 



d 



integrating by parts. 



d_ 

dx \ip 



dx + 



■0 dx 



d 



dx yip 



dx + lim 



(6.4) 



The integral terms in ()6.4p are non negative and equal zero if and only (p = ip. Let us focus on the 
contribution of the boundary, that consists of two terms. Using boundary conditions, the first term 
gives 

lh'n[<P'^Y_^ = -7|0'(O+)P (6.5) 
For the second term we immediately get 

_ v^'(o-)v^(o+)|</.(o-)|2 -V''(o+)0(o-)|</.(o+)|2 



lim 

£^0 



^(0+)V;(0-) 

Summing (16. Sp and (16. 6p . and using the matching condition for both ^ and <j) we finally get 

V^2(O+)|0(O-)|2 + V^2(o_)|^(o+)|2 _ 2V.(0+)V'(0-)Re((/)(0+)^(M) 



(6.6) 



- lim 

£-S>0 



70(O+)V'(O-) 



|^(O+)<^(O_)_^(O_)0(O+)|- 
7V'(0+)V'(0-) 



Due to the hypothesis ^p{0+)^p{0—) < 0, verified by all ground states, we conclude that the boundary 
term in (|6.4p is non negative, and this completes the proof. □ 



Remark 6.2. An analogous proposition holds in the case, treated in [17j . |16j . |28j of a 5 interaction 
with strength a (no matter if attractive or repulsive). One has that (the meaning of the symbols is 
the obvious one) 

a) Ker '"^ = Span {ip} 

b) L°'" > 

In this case, the boundary term in (j6.4p vanishes. 

Now we prove the spectral features of interest for the operator Lj'^ defined in ()6.2p . 
Consider first the case of the symmetric stationary state ip = ■(/'w'"^'^- We define 



w(^ + 1) 



^'^^"^ cosh.'^{fj,y/uj{\x\ + y)) 

and 

'^2,sym - -7 ■ " ^Osh' {fl^{\x\ + y)) 

where tanh(/i-y/tJy) = 7^^- 



Li:^^ = H, + u ,0, , (6.8) 
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Proposition 6.3. Fixed /U > 0, the operator L^'^^^ has: 

a) A trivial kernel and one simple negative eigenvalue, if < u* ; 

b) A one- dimensional kernel, spanned by the function 

sinh(^V^(|x| +y)) 



cosh {fl^/uJ{\x\ ->r y)) 

where y has been defined in (|5.12|) . and one simple negative eigenvalue, if uj = uj* ; 

c) A trivial kernel and two simple negative eigenvalues or a double negative eigenvalue, if oj > uj* . 

Proof. For shorthand, in this proof we will denote the operators ^I'^m -^2 sym -^i ^^^'^ -^2, 
respectively. Furthermore, the function tpJi'^'^ will be denoted by Tp. 

Consider first the case u < u* . By stationarity of ip, the following identity must hold up to 
higher order terms in w: 

Sui{ip + w) = Sui{ip) + ^iu,Liu) +'^{v,L2v), (6.9) 
for any w = u + iv in Q, with u and v real. 

By Proposition 16. H the operator L2 is non-negative, and by Weyl's theorem on the stability of the 
essential spectrum (see for example [23] or Theorem XIII. 4 in [3lj), one has (Tess(Li) = o-ess{H^+^) = 
[^,+00). Furthermore, Theorem 15.31 guarantees that ^p minimizes the functional S^j on the Nehari 
manifold. Thus, since the Nehari manifold has codimension one, the operator Li has at most one 
negative eigenvalue. 
On the other hand, 

{^P,L,i;) = -IfiXUf/^Xl < 0, (6.10) 

so we can conclude that for u < uj* the operator Li has exactly one negative eigenvalue. 

Concerning the kernel of Li, we recall that the only square-integrable solution to the linear 
differential equation 

a ;(/x + l)(2/i + 1) 

is given, up to a factor, by 



e--+u;c- "^^T,7^T^ = (6.11) 



e(^) = — 1 ' ■ (6.12) 

cosh {^y/ujx) 

Furthermore, there cannot be a solution ^ Span(,^) to equation (j6.1ip such that \C,{x)\^dx < 00 
for some finite a, otherwise, by invariance under reflection of (j6.1ip . the function C{~x) would be a 
solution to (|6.12p too, satisfying \C,{x)\^dx < 00, so we would obtain three linearly independent 
solutions to (16. lip . As a consequence, the possible solutions to the equation 

L^c ^ ^.c ^(^+l)(2/i + l) , ^ Q ^g^g^ 
cosh^(^V^ (I • I + y)) 

are given by S,a{x) = x+ix)Cix + 2/) + '^X-{x)(,{x — y), with a G C, provided that they fulfil the 
matching condition of the 5' interaction. Such conditions prescribe the identity of the left and the 
right derivative at zero, namely ^'(y) = a^'{—y). Since ^' is even, this implies either a = 1 or ^'{y) = 0. 
In the first case, imposing the boundary condition ^i(0+) — Ci(0— ) = — 7Ci(0+) leads to the equation 
2 sinh(/x-y/a;y) cosh.{^^/ujy) = ^/uj[s^uh^ [^^/ujy) — that cannot be solved in y for any /i > 0. In the 
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second case, one has ^^(0) = 0) which is fulfilled, as ^^(0+) = \fZj cosh ^ (^■^/6Jy)[;U — sinh^(^A/a;y)], 



if and only if sinh (/iy^y) = /i. This is equivalent to tanh(^-^/a;y) = y that, owing to definition 

of y in (|6.7p . is verified only for a; = a;*. Furthermore, since zero is a stationary point, must be 
continuous at zero, so a = — 1. 

Thus we proved points a) and 6), and the case cj < is exhausted. 

In order to prove point c), let us write the spectrum of L\ as 

(t(Li) := {i/i, . . . U {n, . . . ,Tm} U [w,+oo), 

where Vi < Uj < and > > for any i < j and any / > m. The boundedness and fast decay in 
X of the last term in the l.h.s. of (|6.13p ensures that both m and n are finite. Moreover the essential 
spectrum coincides with the one of Li, thanks to the Weyl's theorem again. 

By ()6.10p . we know that n > 0. Denoted by Pa the orthogonal projection in on the eigenspace 
associated to the eigenvalue a, we define the following operators: 

- P- is the projection on the space 0^ P^. ; 

- P+ is the projection on the space Pr^ ; 

- Pc is the projection on the space associated to the essential spectrum [w, +oo) of Li. 

Let us suppose that n = 1. Then, denoted by Fi the quadratic form associated to the operator 
Li, there exists at least a non vanishing combination rj oi ip and that satisfies Fi{r]) > 0. Indeed, 
denoted by V'l the only (up to a phase) normalized eigenfunction associated to the eigenvalue —I'l, 
we define the function 

(V'l^C-i) , , . 

V ■= — n 7V^ + ^-i- 

(V'l,^') 

Notice first that, since by (j6.10p Fi^ip) < 0, ^ and ipi cannot be orthogonal, therefore rj is well 
defined. Furthermore, ip and are linearly independent as they have different parity, so it must 
he 7] ^ 0. Since {ipi,r]) = 0, tj has no components in the negative part of the spectrum of Li, so 
Fiiv) > 0. 

But this is not the case. Indeed, for a generic combination (p = a'lp + /3^_ , we have 

Fi(0) = |a|2Fi(V^) + |/3|2Fi(C_i) + 2Rea/3(LiV,e-i> = |a|2Fi(V') + |/3pF7'"(e-i), (6.14) 

as the mixed term vanishes, being the scalar product of an even and an odd function. 

Now we compute Fi(^_i). We notice that, due to the continuity of the term related to the point 

interaction vanishes, so, after integrating by parts, we get 

+{_i(o)({'_,(o-)-!;'_i(o+))] 

= i-mi'-io-) - Ciio+)), 

where we used the fact that, by definition of 

-ei(x) + .C_i(x) - "gy+^^/,^ + ^\. e^i(x) =0, Vx / 0. 
cosh {iiy/uj{\x\ + y)) 
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Then, one can directly compute 



-ill 

7 



7^0; 



which is negative if and only if a; > w*. Thus, as a consequence of (|6.14p . beyond the bifurcation 
frequency oj* , for any linear combination (p of and ip we have Fi{(p) < 0, that contradicts the 
hypothesis of having only one simple negative eigenvalue in the spectrum of Li. 

In order to prove that actually either n = 2 or n = 1 and vi is a double eigenvalue, we prove 
that minimizes S^^ on the Nehari manifold with the additional constraint 92(0+) = —ip{0—). To 
this aim, we first observe that, if 99 E Q fulfils 92(0+) = —ip{0—), then 



1 II / ii9 II ii9 

2"^" 2"^" 2^+2 



7 



I 'Il2 , II ||2 \|| ||2u+2 

\ip II + ijo\\ip\\ — AII99" 



.16) 



2^+2--iv^(o+)r 



Consider the unitary transformation C/j of the space Q, defined by 

ip^{x) := {U^f){x) := e{x)(p{x) 

and notice that, if 9?(0+) = —ip{0—), then ^p^ belongs to -ff^(M), so the minimization problem is 
mapped into the problem of minimizing the functional 



1 



/ ||2 



UJ 

+ 2 



A 



2iirjii 2"^«" 2/i + 2 
among the functions in (R) that satisfy the constraint 



(6.17) 



A 



0. 



Problem (|6.17p corresponds to the issue of finding the ground state for a nonlinear Schrodinger 
equation in the presence of a 5-type defect of strength —4/7. By [17J and [H], we know that the 
solution reads 



so we obtain, still up to a phase. 



A 



1 

2m _1 

cosh 



\x\ +y) 



and, as a consequence, ip minimizes the action on the Nehari manifold with the additional condition 

V'(o+) = -V'(o-). 

It remains to prove that such a constraint has codimension two. We denote the constraint by 



(6.18) 



By the operator U» such a constraint is mapped to 



Mi := {ipeH\M.), /^,s((^) = 0}. 



.19) 
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It is well-known ([17j. [16| ) that A^j has codimension one as a subspace of i?^(R). On the other 
hand, since any function C in Q can be decomposed as 

C = ^(C(0+)-C(0-))6(-)e-H+C, 

with C £ H^O^), we have 

Q = i7i(M)eSpan{e(-)e-|-l}, 

so H^(W) has codimension one as a subspace of Q. Therefore, has codimension two as a subspace 
of Q. Thus, by unitarity of C/jj, has codimension two as a subspace of Q too. 

We then proved that the negative space of the operator Li has dimension at most two and that 
for w > a;* it equals exactly two. The proof is concluded. □ 

Now we consider the case of an asymmetric ground state ipJ!}' v^fi ^ define 

,7- - U ^ w(/i+l)(2/i + l) 



cosh^ (/i^ {x + x+{x)y2 - X- {x)yi)) 
and 

-7,0; _ TJ , . 



^ cosh? {x + x+{x)y2 - X-{x)yi)) 

where tanh(^-y/a;yj) = tj, and ti,t2 are the unique positive solutions to (15.20 with ti < t2- 

Proposition 6.4. Fixed /i > 0, for any uj > u* the operator Ll'^^^ has trivial kernel and one 
simple negative eigenvalue. 

Proof. For shorthand, in this proof we will denote the operator i^'asym -^i- Furthermore, the 
function will be denoted by tp. 

By Theorem 15.31 if w > w* then ip is a local minimizer for the functional S^j on the Nehari 
manifold. As a consequence, one can prove that Li has one simple negative eigenvalue by following 
the proof of Proposition 16.31 up to (j6.10p . 

Concerning the kernel of Li, one can follow the reasoning carried out in the proof of 16. 31 through 
(j6.1ip . (|6.12p and conclude that the only solutions to the equation Li^ = can be given by (,aix) = 
X+(,{x + 7/2) + cix~(,{x — yi), where o is a complex number, provided that fulfils the matching 
conditions at zero, that translate into the system 



2+1 ~ ^ 2+1 

cosh (M\/iiJ2/2) cosh {^^/Zjyi) 

ain\^{^l^y2) . sinh(^v/^yi) _ &\nY^{^ly^y2) 

1+1 " 1+1 ~ — Tv"^ 5^"T 

cosh cosh ''(/i-y^S/l) cosh ^(/Jv^J/2) 



(6.20) 



Expliciting a from the first equation, plugging it into the second, and denoting as customary ti 
tanh(/iy^yj), we get the equation 

^ 7 — r-rTT2 = -7vw, 



that, using the second equation in (j5.2p . gives 



1 - i-tj 



tiifi - {fi + l)ti) tiifi - ifi + l)ti 
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Finally, by the first equation in (|5.2p one gets 



Such a problem translates to the problem of finding ti and t2 such that < ti < i < t2 < 1, and 

= 9Ht2), where ^(t) := ^^'(0 and f{t) = t^^ - t^M+s. 
The problem was treated in the last form in the proof of Theorem 15.31 from formula (|5.19p up to 
formula (I5.26p . The conclusion is that it has no solutions, so none among the functions lies in the 
kernel of Li, which is therefore trivial. This concludes the proof. □ 

6.2 The sign of 

Proposition 6.5. Given fJ- > 0, the sign of the second derivative of the function d, defined in ()2.13p . 
is determined as follows: 

1. ifO < /i < 2, then d"{uj) > for any co G {ujo,uj*) U (a;*,+oo). Furthermore, < d"{cL)* + 0) < 
d"{uj* - 0); 

2. if fJ.> 2, then d"{u]) > for u; £ (wq, oj*). Besides, d"{uj* - 0) > 0; 

3. there exists /u* G (2, 2.5) such that 

(a) If pi < fJ-*, then d"{oj* + 0) > 0, so there exists f^i(^) > w* such that d"{oj) > for any 

(h) if fi = fi*, then d"{uj* + 0) = 0; 
(c) if fi> fi*, then d"{uj* + 0) < 0; 

4- if fJ- > 2, then there exists 0J2{pi) > w* such that, if u > W2(^), then d"{uj) < 0. 

Proof. First, we notice that, given n > 0, u > ujq, and denoted by ip^^ a solution to the problem 12.31 
corresponding to the chosen value of /x and to, one has 

Indeed, from definition (j2.13p . using the stationarity of ipuj we get 



Expliciting M(ipi^) one obtains 



d'iio) 



\ J 2^ 



1 1 rl 1 

{l-t^)^-^dt+ / {l-t^)^~^dt 



(6.21) 



provided that the change of variable t = ianh.{iiy/oj{x — xi)) {t = tanh(//\/tJ(a; — X2))) has been 
performed in the first (second) integral. The lower bounds in the intervals of integration are defined 

by 

CM = { , ^ = 1,2 (6.22) 

ti, UJ G {ijj* , +00) 
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where the couple ti,t2 is the unique solution to the system (j5.2p such that ti < t < t2, where 
t = 



Differentiating ()6.2ip yields 
d"{oj) = 



1 1_ 3 

/U + IXmo;'' 2 1/1 1 



A 



2fi \\n 2 



1 1 /■! 1 



-OJ 



Let us denote 



(I) := 



(11) :- 



C((a;)(l - C?H)^'' + C^(^)(l - ClH)^~' 



3.23) 



(6.24) 



1. < n<2. 

The quantity (/) is positive for any u > ujq. Moreover, by (j6.22p . for uj G {ujq,uj*) the quantity {II) 
can be explicitly evaluated as 



(11) 



.25) 



which is positive too. Therefore, by (j6.23p and (|6.24p . we conclude that d"{uj) > 0, if t^o < < 

To determine the sign of (//) for a; > w*, it is convenient to distinguish between the cases /U < 2 and 

fj, > 2. Let us start with fi < 2. Rewriting the first equation in (15. 2p as 



and differentiating with respect to cu, we obtain 



2\ --1 



fl 



where we neglected in the notation the dependence on u. Therefore 



(11) 



We prove that such a quantity is positive for fi < 2. As ti < t < t2, this reduces to prove 

{fi + l)4-f^t2 > fih - {fi + 1)4. (6.26) 
To this aim, we define the function 

m := [(^ + I)t3_^t]2 = 1^4-4M[^/(^)]2^ 

where / is defined as in (j5.14p . namely by f{t) := — By the fundamental theorem of the 

integral calculus, for any t G [0, 1] 



ii-fi)s'-^^[f'{s)f + ^s^-'^ns)ns) 



ds, 



(6.27) 
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where t = yj^- Now we proceed like in the proof of theorem 15.31 formulas (|5.20p - (j5.25p . First, we 

define the function si as the inverse of / in the interval [0,t|, as well as the function S2 that is the 
inverse of / in [i, 1]. Then, performing the change of variable u = /(s), ()6.27p gives 

T{t) = / Y.{si{u))du (6.28) 



where i = 1 if t G [0, t], i = 2 if i G (t, 1], m := /(f) = ^^^^^y+i , and 

S(s) := (1 - fi)s'-'^ns) + ^.^-^'^r (.) = f.s'-''' - 3(^ + l).^-2^ 

Consider the case h £ [t/^/3,t). The study of the sign of S and S' shows that S is negative and 
strictly decreasing in (t/\/3, 1), for any u G [f{ti),m) one has 

> S(si(n)) > m = -^J^^J^ > S(s2(n)), 
and therefore, denoting a = f{ti) = f{t2), 

T{ti) = - S(si(u))dn < - / S(s2(n))dn = T{t2). (6.29) 



Second, consider the case ti < t/Vs. Write r(ti) as 

rfit/V3)) rfitj) 
r(ti) = / Y.{si{u))du+ / J:{si{u))du. 

Notice that the first integral in the r.h.s. is positive, while the second is negative, owing to the facts 
that S is positive in {0,t/^/3) and that f{ti) < f{i/\/3). Then, denoting by t2 the only point in 
(ti, 1] such that f{i/V3) = /{h), by (lOOl) . 

r{ti) <r{t/V3) < Tih). (6.30) 

Furthermore, since ti < i/y/S, it must be ^2 < ^2) and since S is negative in the interval {i2,t2), 

we obtain r(t2) = - Z/^*^-* T,{s2{u))du + r(t2) > r(t2) that, together with (fHIHOj) and ([6:291) yields 

r(ti) < r(t2) for any ti G [0, t), which is equivalent to (|6.26p . So we proved {II) > if /i < 2, that, 

together with the positivity of (/), proves that d"{uj) > for any a; > wq, w / a;*. 

In order to complete the proof of point 1., we need to evaluate d"{uj* it 0), and to compare them. 

By (f6:23]l and (fHIM]) . this amounts to comparG tli6 value of the two limits lini^^^^* ±o[{I) + iII)]. From 

(|6.24p it is clear that (!) is continuous at lo* , where it takes the value ( f — 1) .f^r~P~i^ —t^)i^ ^ dt, 

^ V 

so we reduce to study the limits of the term (//) only. The left limit is immediately given by (|6.24p 
and ()6.25p . and reads 

hm (11) = > 0. (6.31) 

Computing the right limit is more complicated. Differentiating both equations in ()5.2p one can 
express the couple (^1,^2) as a function of the couple (ti,t2)) namely 

^ int2)44 ^ if'{h)tjti 

2^{tir{t,) + tir{t2)y 2v^(tfm) + tim))' ^ ^ ^ 
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where we neglected in notation the dependence of ti and t2 on the variable uj and denoted f{t) = 
^2fi _ ^2^+2 _ Yrom (f02D . ^22i), dS^D, and since, by ([O]), f{ti) = f{t2), one gets 

^''^ - ^ ^''^ tim)+tint2) ■ ^'-''^ 

Here, the only non trivial factor is the last one, in which both numerator and denominator vanish 
as oj goes to oj* from the right. In order to compute such a limit we pass from the variable uj to the 
variable ti. In other words, we consider t2 as a function of ti. 
We define the functions 

N{tr) := tl-'^nh) + tl~'^{h)f'{t2{h)) 
D{h) := tlf{t,) + 4{t,)f'{t2{h)) 

and provide a Taylor expansion near ti = t = -^^j^ for both of them. One immediately gets 

iV'(ti) := (2-2^) [t}-''^/'(ti) + 4"'^/(t2)i2] +t?""^r(ti) + t2'''^/"(t2)i2 
Z)'(ti) := 2ti/'(ti) + t2/"(t,) + 2t2/'(t2)i2 + i^/"(i2)i2 

where we used the notation {2 := ^(ti)- To evaluate A^'(t) and D' {t) we must then compute ^(f)- 
From /(ti) = /(t2) it follows 

^■2 = 

By de I'Hopital's Theorem 



H ^ (0.35) 



lim ^2 = liin ^ 



ti^t-O ti^t-O f"{t2)t2 

from which one immediately has (lim^j_!.j ^2)'^ = 1- Now, since /'(ti) > and f'{t2) < 0, it must be 

lim t2 = —1. 

As a consequence, A^'(f) = D'{t) = 0. Further differentiating and D, and recalling that f'{t) = 0, 
one obtains 

N"{t) := 8(l-/u)t^-2/^/"(t)+2?-2M/"(t)+t2-2Mj//(i-^^(i-^ 

•,r2j:////T\ I +2i-///7\" ''2 



:= 8t7"(t) + 2t^r(t)+t'/"(t)3lf©- 



By ([05]) 



^ ^ r(ti)(/o^(t2(ti))-(/o^(ti)r(t2(ti)) 

so, again using de rHopital's theorem, 

t,^-o 3f"(t} ' 
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It then follows 



N"{t) := 8{1- ii)f-^f'f"{t) + fit) 

3 

D"{t) := m"{t) + U^fit). 



So, going back to (|6.33p . we get 



lim {II) 



12/^1'^ fl-\i)^V: - ^)i'-'^f"{i) + P-^^f"'{t) 



Qtf"{t) + Pf"'{t) 



which, recalling that t = \/ j;^ and the definition of the function /, gives 

lim (11) = > 0. 

a;->w*+o _^ 1)71-2 4// + 5 



(6.36) 



Comparing (j6.3ip and (j6.36p . and observing that the existence of such limits implies the existence 
of the left and right derivative and gives their values, one completes the proof of point 1. 

2. fi > 2, bjQ < oj < uj* . 

In this case, term (/) is negative, so one must compare its size to the size of (//). 
Prom ([6T22|) and ([03]) we know that 



where 



Then, 



r(u;) 



d"{uj) 
2-fi 



2u 2 

7v^ 



{l-t^)T^~^dt+ ^ 



r{u}), 



r'ioj) 



3 

JUJ2 



1-- 1 ( 1 



7VW 



4 \ M 



7^0; 



< 0. 



We estimate the first term in the r.h.s. of (j6.38p . for u = lo* , as 



> 



2-fi 
2/x 



:i-t^)^''dt > 



2-^1- 



M+1 



2- ^l 



M + 1 



2^ 



1 + 



2(1 + a/^)(/^ + 1) 



(6.37) 



.38) 



So, 



r{uj*) > 



1 



2{^l + l) 



2 - /i 



1 + ^/^1 



+ V/x(/i + 1) 



2(/i+l)M 



Thus, since r is monotonically decreasing, we have r > for any r G {u}o,u}*]. 

Finally, since liuica^uj* -0 d" (u) exists and can be recovered by putting uj = uj* in (I6.37p . one 
obtains that d"{uj* — 0) equals such limit and so point 2. is proven. 
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5. ^ > 2. 

In order to prove points 3 (a,b,c) we need to evaluate d"{uj* + 0). We compute liuiui^uj* +o d" (co) . 
Then, 

lim d"{(jj) = lim (/) + lim 
with (/) and (//) defined in (j6.24p . and, by a direct computation, 

hm (!) = (--l] I {l-t^)^'^dt < 0. 



Exactly as in the proof of point i., one finds that the right limit of the term {II) is given by formula 
Now, according to (|6.23p and (|6.24|) . the sign of d"{uj* + 0), is given by the sign of the function 



M/.) := hm [(/) + (//)] = (--l) [a-t'y-'dt 



(/i + l)M-2 4^ + 5' 



Such a sign is obviously negative for /i > |, so we restrict to 2 < /i < |. Then, 



(it 



M+l 

1 1 



^2V^(// + l)i+^4// + 5 ;j "^(^ + 1)1-^ 4// + 5 

^(// + 1)^-M 



- 30 



(4/i + 5)2 

(6.39) 



For 2 < // < I the only non negative term in (j6.39p is the third one. By elementary computation we 
find that the sign of the sum of such term with the last one is negative if and only if 

-16/i^ + 12/4^ + 25 < eO/i^. 

Since for // = 2 the inequality is verified, it must be verified for any // > 2, owing to the fact, 
easy to check, that the l.h.s is a monotonically decreasing function, while the r.h.s is monotonically 
increasing. As a consequence, w is monotonically decreasing for // E [2, |], so there exists a unique 
value of // that makes d"{u}* + 0) vanish. Denoting it by //*, we complete the proof of points 3 (b,c). 
Point 3 (a) follows by continuity from the fact that, if /i < /i*, then d"{uj* + 0) > 0. 

4. /i > 2, w — )■ 00. 

Point 4. follows from the asymptotics of ti and t2 as uj goes to cxd, in the region t2> ti > 0: 

h = -^ + oiuj'h t2 = 1 ^ Voioj'^) (6.40) 

A = ^ + o(a;-i) t'2 = , 2a u+i +^i^''-')- (6-41) 

27a; 2 27^^^^*+^ 
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Let us prove such asymptotics. The condition t2 > ti selects, among the solutions to the system 
(j5.2p . those belonging to the set T2 defined in (|5.16p . It is immediately seen that, in such region, 
ti — )• and ^2 — ^ 1 as w — )• 00. As a consequence, from the second equation in (|5.2p one gets 

lim 'J^/u}tl = lim — "^^^ ^ — = 1, 

SO the first formula in ()6.40p is proven. 
From the first equation in ()5.2p we have 

lim tV^f^il-tl) = lim t2 ^''(1 - *i) = 1, 

and by the first of (tO0]l 

1-4 = -,^ + o{u>-^') 

and the second identity in (|6.4Up immediately follows. 

In order to prove (j6.4ip we differentiate both equations in (j5.2p . obtaining 

fl f-"^ _L +'+-2 _ 2_ 

Expressing t'2 from the second equation and plugging it into the first one we get 



(6.42) 



2tf +^(^ - (/i + l)tl) + 2tf +^(/. - (/i + l)t\ 



that converges to —2!^ as a; goes to infinity, and so the first formula in (j6.4ip is proved. To prove the 
second one, it is sufficient to use the first equation in (I6.42P and by (|6.40p . (j4.9p . the first in (j6.41jl 
one finally obtains 

lim u^'+^t' = -At 

!>00 27 

and the proof is complete. □ 

6.3 Stability and instability of the ground states. Pitchfork bifurcation 

We recall the definition of orbital neighbourhood and orbital stability. 
Definition 6.6. The set 

U^{^) := G Q, s.t. inf ||^ - e'Vllq < v} 

0e[O,27r) 

is called the orbital spherical neighbourhood with radius r] of the function (p. 

Definition 6.7. We call orbitally stable (in the future) any stationary state (j) such that for any 
e > there exists 5 > such that 

inf IIV^ - e^VllQ < 5 ^ sup inf \\m - e'^h < e, 

6'e[0,27r) t>0 De[0,27r) 

where 'ip{t) is the solution to the problem (13. ip with ip as initial datum. 
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Figure 2: Hypocritical and critical case (i.e. /U < 2). All ground states are stable. 



Definition 6.8. We call orbitally unstable any stationary state that is not orhitally stable. 

Proposition 6.9 (Stability and instability of ground states). Consider the ground states of the dy- 
namics described by (jl.ip . defined as the solutions to problem \2.4\ and explicitly computed in Theorem 
nOl Then, 

1. If < fj, <2, then for any u G (ojq, +oo), oj ^ lo* , all ground states are stable. 

2. If n> 2, then 

(a) there exists uji > u* such that all the ground states with frequency u < loi are orbitally 
stable. 

(b) There exists uj2 > uji such that if u > UJ2 then all ground states with frequency to > U2 are 
orbitally unstable. 

Proof. Points 1 and 2 (a) follow from Theorem 2 in [21]. Indeed, notice that Assumption 1 in such 
theorem is proven by Propositions 13.3 1 and 13.41 while Assumption 2 follows from Proposition 15.11 and 
Theorem 15.31 Furthermore, owing to propositions 16. li 16.31 a), b), and 16.41 Assumption 3 is verified 
for all ground states. Therefore, in order to establish orbital stability, it is sufficient to remark that, 
for the considered cases. Proposition 16.51 establishes d"{Lo) > 0. 

Case 2 b) follows from Theorem 4.7 in [21j, as we know from Proposition 16.51 that d"(uj) < 0. 
Thus, the theorem is proven. □ 

Remark 6.10. By Theorem 3 in pT] we have that, under the hypotheses of Theorem [621 the ground 
states solve Problem 12.21 they minimize at least locally the energv [2181 among the functions with 
the same L^-norm. 

Theorem 6.11 (Pitchfork bifurcation). Given fJ. > 0, if u > uj* then the stationary solutions ijju ^'^ 
defined in Theorem 15.31 are orbitally unstable. 
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Figure 3: First hypercritical case (i.e. 2 < fx < ix* < 2.5). Symmetric ground states are 
stable. Immediately after bifurcation, the two newborn asymmetric states are still stable. At large 
frequencies, they become unstable. 



t2 



W2 



U) > U) 



W2l 



Figure 4: Second hypercritical case (i.e. // > Symmetric ground states are stable. Change 
of stability occurs immediately after bifurcation. 
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Proof. Closely mimicking the computation that led to formula (|6.23p . one gets 



As uj > Lo* , we know by proposition 16.31 we know that the number of negative eigenvalues of the 

operator L'^''^ equals two. 
from Theorem 6.2 in 1221. 



operator L'^''^ equals two. If — "^^^^ i^) is positive (for instance if ^ < 2), then the result follows 



On the other hand, if — ^^^^ i^) is negative, then we restrict the problem to the space Qa 

of antisymmetric functions lying in Q. First, from the explicit knowledge of the propagator of the 
Schrodinger equation with a 6' interaction, represented by the integral kernel (see [1]) 



:{\x\ + \y\-u)'^ 



g * T (x v) = h €(xv) I — - — - / e I — -e ^ 

27 7o 7 

it appears that, denoting g{x) = g{—x), one has 

Let us consider the problem (jl.ip . and initial data tpQ such that ipQ = — V'o- Then, applying Duhamel's 
formula to (jl.ip . one finds 

Jo 

so that 

Jo Jo 

It follows that ipt solves (jl.ip with —ipo as initial data. Since (jl.ip is invariant under multiplication 
by a phase factor, and since the solution is unique, it must be tpt = —tpt, and so we have that (jl.ip 
preserves the antisymmetry, and therefore the evolution problem (jl.ip is well-defined in Qa- 

Now, as already remarked in the proof of Proposition [631 the functions are the minimizers 

of among the antisymmetric functions belonging to the Nehari manifold. The time-evolution 
operator, linearized around them, admits a simple negative eigenvalue, so, by Theorem 4.1 in |21j . 
the stationary state V'S'"^' is orbitally unstable in Qa, so, a fortiori, it is orbitally unstable in Q. 

To complete the proof, we recall that the instability of ip^~^' in the case — '^^^ji {'^) = can 

be established by the argument in [12]. □ 

We are then in the presence of a pitchfork bifurcation, that can be depicted as in Figure 5. 

We leave open the issue of optimizing oji and 0J2, and possibly getting uji = uj2- In other words, 
we do not know whether, for a frequency beyond the bifurcation, each ground state undergoes only 
a change of variable or more. 

Finally, we don't treat here the problem of determining the stability features of the ground states 
at the bifurcation frequency. This will be the subject of a forthcoming note. 

Acknowledgements. The authors thank Andrea Sacchetti for many discussions during the 
initial step of the project, that helped us to focus the problem. R. A. is grateful to Cinzia Casagrande 
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Figure 5: Pitchfork bifurcation. The branch of the (anti-)symmetric ground states can be con- 
tinued beyond bifurcation, for any value of oj, but the corresponding stationary states are orbitally 
unstable. By the same modification, i.e., adding the unstable straight branch from the origin up to 
to Figures 3 and 4, one obtains the pictures that correspond to the cases 2 < /j, < n* and 



for some explanations on algebraic curves and to Sergio Conti for an illuminating discussion on the 

properties of the symmetric stationary states beyond bifurcation. 
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Methods for Nonlinear Differential Equations. 
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